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Abstract 

The statistical properties of failure are studied in a fiber bundle model 
with thermal noise. We show that the macroscopic failure is produced by 
a thermal activation of microcracks. Most importantly the effective tem- 
perature of the system is amplified by the spatial disorder (heterogeneity) 
of the fiber bundle. The case of a time dependent force and the validity of 
the Kaiser effects are also discussed. These results can give more insight 
to the recent experimental observations on thermally activated crack and 
can be useful to study the failure of electrical networks. 

PACS: 05.70.Ln, 62.20.Mk,61.43.-j 

1 Introduction 

Material failure is a widely studied phenomenon not only for the very important 
technological applications but also for the fundamental statistical aspects, which 
are not yet very well understood. Many models have been proposed to give 
more insight into the statistical analysis of failure. Among the most studied 
ones we can mention the fuse and the non-linear spring networks jl], ^, ||, |], [|, 
[| 0], which can reproduce several features of crack precursors experimentally 
observed in heterogeneous materials subjected to a quasi-statically increasing 
stress H H [it], [n], [l^, [l3| . Specifically the power law behaviour of the acoustic 
emission observed in several experiments close to the failure point. However 
the non- linear spring networks and the other related models, in their standard 
formulation, are unable to describe the behaviour of a material subjected to a 
creep-test, which consists in submitting a sample to a constant stress till the 
failure time. Creep-tests are widely used by engineers in order to estimate the 
sample lifetime as a function of the applied stress. Modified fuse networks have 
been proposed to explain the finite lifetime of a sample subjected to a constant 
stress. In ref. [El for example, the Joule effect has been considered to explain 



delayed crack. Instead in ref.|| [| phenomenological force dependence of the 
fiber lifetime has been proposed, but not very well justified from a physical point 
of view. These models do not explain the recent experimental results on 2-D 
crystals [O, gels|ll| and heterogeneous materials Q. 

These experiments show that the lifetime r of a sample, subjected to an 
imposed stress P is well predicted by the equation 

/ T d Y (d-i) \ 
t = t exp a-— -^-—r (1) 



K T e ff P(2rf-2) 

where r Q is a constant, T the surface energy, Y the Young modulus, K the 
Boltzmann constant a a constant, which depend on the geometry, T e ff is an 
effective temperature and d the dimensionality of the system. Eq.|l| was first 
derived by Pomeau and Golubovic jlT], for d = 2 and d — 3. It has been 
generalized to any d in ref. |l5| . The main physical hypothesis behind eq.[l] is 
that the macroscopic failure of a material is produced by a thermal activation of 
micro cracks [|~7j [lg]. In the original Pomeau's theory JlTj T e ff of eq.0 coincides 
with the thermodynamic temperature T while, experimentally, T e fj >> T jl4], 
|l5[ [ll| . Eq. ([!]) is based on the idea that the sample failure is due to the thermal 
nucleation of one defect (microcrack) initially present inside the material fl7f. 
However, experimentally it has been observed that the macroscopic failure is in 
some way related to the nucleation of many defects; this important point will 
be discussed at the end of this paper. Using reasonable values for the material 
constants in eq.[j] jl6| one sees that thermal fluctuations are too small to activate 
the nucleation of microcraks in times r measured in the experiments. It has been 
measured that the temperature needed to have the measured lifetimes r should 
be of order of several thousands Kelvin. Specifically it has been estimated that 
for wood T eff > 3000K Jll , in 2D crystals lOOOif < T eff < 2500K @ and 
T e ff > 10 10 K for gels |15[. In several materials the lifetime r of the sample 
is not influenced, in the limit of experimental errors (20 %), by a variation of 
the temperature T from 20 to 90 °C. In contrast, experiments on 2D-crystals 
[ pT[ show that t depends on T. Furthermore the lifetime statistics is close to 
Gaussian in ref.@, @. Thus eq.|l| seems to give the right dependence on P 
of the sample lifetime but there are a lot of experimental observations that are 
not completely in agreement with a standard nucleation problem. To explain 
these results it has been supposed Q that the lifetime r of the sample depends 
on the heterogeneity of the material and that disorder in some way enhances Q 
thermal fluctuations so that the nucleation time of defects becomes of the order 
of the measured ones. 

Motivated by these experimental observations and to check the validity of 



our hypothesis we proposed 19, 20(] a model based on a democratic fiber bundle 
model (DFBM) with thermal noise. The DFBM is the simplest spring network 
proposed a long time ago by Pierce [|2lf to study cable failure. The DFBM 



has been studied by several authors 111, El tj, E2L E3|, and it turns out to be a 



1 The term "enhances" has to be understood in the sense that the effect of the disorder is 
such that the system behaves like if the real temperature is bigger 
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quite realistic model for long flexible cables, low- twist yarn and more gener- 
ally for composite heterogeneous materials. In a classical DFBM, which is a 
deterministic model, it is not possible to observe stochastic phenomena such 
as nucleation. Introducing thermal noise into the system we transformed the 
DFBM in a stochastic model, so that processes such as nucleation become ob- 
servable. 

The purpose of this paper is to extend the results of our previous work p9l ^Oj 
and those obtained analytically by Roux pi[| . We will show that, adding to the 
spring network a noise, which plays the role of a temperature, it is possible to 
reproduce the behavior of a material subjected to a creep-test. Furthermore 
we can prove that, for such a process, the noise is "amplified" by the network 
disorder. 

The paper is organized as follows. In section 2 the modified DFBM model 
is introduced. In section 3 the statistical features of the model are studied 
analytically. Specifically we compute how failure time r is related to the disorder 
of the material, when the bundle is subjected to a constant force. In section 4 the 
results of section 3 are verified numerically. In section 5 we study the lifetime in 
the general case of a bundle subjected to a time dependent force F(t). We also 
propose a formula, which allows us to estimate the failure time r by knowing 
the whole history of the medium, that is to say, the time dependence of the 
imposed force F(t). In section 6 we give a comprehensive explication of some 
facts related to Kaiser effect, with a comparison between the experimental data 
and the numerical results obtained with our model. Discussion and conclusions 
follow in section 7. 



2 The model 

We simulate the behavior of an elastic heterogeneous medium subjected to an 
imposed tensile force F using a DFBM with thermal noise. For this purpose, we 
model the medium as a system of N parallel clastic fibers whose extremities are 
fixed on a rigid support as shown in fig. [l| This model, widely studied in the 
quasi static regime [jj], ||, f|, is equivalent to N springs in parallel subjected 
to a total tensile force F. Specifically we have studied the model by using the 
following rules: 

I The external applied force F produces a local force fi on each fiber. F is 
democratically distributed in the net: F — Yli=i fi- 
ll The local force fi on the i-th fiber produces a local strain e,. Being the 
media elastic, force and displeacement are linked by the Hooke's law: 

fi = Y ■ e t (2) 

where Y is the stifness, which is assumed to be the same for all the fibers: 
Yi=Y, 
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III The strength of each fiber is characterized by a critical stress / c (i): if at 
time t the i-th fiber local stress fi(t) is greater than a critical stress f c (i) 
the fiber cracks, and the local force falls to zero at time t + 1. Further, 
we assume that in this process some energy ej, proportional to the square 
of local force /i, is released. For sake of simplicity we will assume Ci = 
0.5 ff/Y. The critical stress f c (i) is a realization of a random variable 
that follows a normal distribution of mean < f c > and variance KTd'. 

f c (i)~N d {<f c >,KT d ) (3) 

We call Nd the disorder noise. 

IV Each fiber is subjected to an additive random stress Afi(t), which follows 
a zero mean normal distribution of variance KT: 



Afi(t)~ N T (0,KT) (4) 

being t the time. We call Nt the thermal noise. We assume that Afi(t) 
is a white random process, which is independent in each fiber, i.e. the 
correlation function E [A/j(ii) • Afj(t 2 )] — if t\ ^ t 2 or i ^ j. A time 
step of the model corresponds to the application of the noise to all the 
fibers. H 

The first three items are those used in the standard formulation of the fiber 
bundle model. A new feature, which is similar to a thermal activation process, 
is introduced in [IV] to explain the dependence of the failure time on a con- 
stant applied stress. The model has the following properties. We see from [I] 
that there exists a long-range interaction between the fibers, that is, under the 
specified boundary conditions, the breaking of some fibers produces an increase 
of the local force on the other ones. Indeed if a number n(t) of fibers is broken 
at time t then the local force on each of remaining fibers is: 

/*(*) = IV^+^t) (5) 

where f Q = F/N is the initial force per fiber. Procedure [II] specifies the elastic 
behaviour of fibers. Procedure [III] models the heterogeneity of the medium: 
if there is no disorder noise (i.e. KTd = 0) all the fibers are strictly equal, 
and if there is no thermal activation (i.e. KT = 0) the medium breaks down 
instantaneously when F — Nf( c \ Conversely, if KTd > fibers begin to crack 
before the final break-down of the media. The assumption that all the disorder 
in the model appears in the strength distribution rather than in the elastic 
constants may be argued by noticing that the effective elastic constant of a 
single fiber is essentially the average of the local elastic constant along the fiber, 
while the strength is determined by its weakest point || . A thermal activation 
process is introduced in [IV] to explain the dependence of failure on time; while 
[I-III] are standard procedures, [IV] is new with respect to DFBM. 

In this paper we want to study the failure time r of the sample as a function 
of /, KT and KT d . 
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3 Some analytical results 



3.1 The energy barrier 

Eq.[j] has been obtained by Pomeau by considering that the final fracture of a 
material is produced by a thermally activated nucleation process of microcracks. 
The energy barrier of this processes has been estimated by using the Griffith 
criterion |Q for the stability of a fracture in an elastic medium. In this sec- 
tion we want to construct a similar criterion for the bundle by using the same 
procedure used in [|l7| p6| . More specifically we compute the energy barrier 
that the system has to overcome in order to break the bundle when n fibers are 
broken and a total constant force F is applied. We begin by noticing that the 
strain energy in the bundle with a constant tensile force is Ue = F 2 X/2 where 
A = 1/[Y (N — ri)] is the bundle compliance. The change in the potential energy 
due to the damage extension is AUa = —F 2 A\. Thus the mechanical energy 
is U M = U A + U E = -F 2 /[2 Y (N - n)] + F 2 /{2 Y N). The energy U f used 
to break n fibers can be computed by considering the energy needed by n fibers 
to reach the elongation corresponding to / c , that is Uf = f 2 /(2 Y)n. Thus the 
total energy ( for Y=l) is 



Etot — 2 



F2 F 2 



N N-n 

At n = n c = N — F j f c this energy has a maximum: 

f 2 N ( F 



(0) 



E max - c 2 \^ N ^ lj (7) 

This means that for n > n c the force per bond is larger than f c . Thus the 
bundle is unstable and will break. Following Pomeau's activation model the 
minimum energy per fiber E max /N has to be compared to KT to estimate the 
lifetime of the bundle. Thus the lifetime of the sample has a different functional 
dependence than eq.H: 



r exp 



fc (f 



2KT V/, 



- f 



(8) 



where / = F/N. Comparing eq.|l| with eqj| we see that the functional 
dependence of r on / is different for the bundle (eq.||) and for a solid (eq.|l|). 
This difference is obviously due to the different geometry of the tensile force 
redistribution in the bundle and in the solid. In section 4 we will show, using 
numerical data, that eq.|| well describes the behaviour of the lifetime of the 
bundle, measured in numerical simulations. However eq.^| can be obtained by 
a statistical analysis of the DFBM with noise. A similar analysis has been 
performed by Roux p4| and in the next section we extend his analysis. 
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3.2 Statistical analysis 

Wc describe here only the main results. The details of the calculation can be 
found in the appendix. 

3.2.1 The lifetime of the homogeneous fiber bundle 

We first consider the homogeneous case where the threshold force f c = 1 is the 
same for all of the fibers. Let P(rj) = ^ 2 ^ KT CX P ( — ■n?TJ ^ x * ne probability 
that the thermal noise has an amplitude larger than r\. The mean force present 
on each fiber is / = fojfz^- 

It can be proved that in the limit N — > oo the ratio becomes independent 
of TV (see appendix). Thus we can introduce the function <f>(t) (fraction of broken 
fibers at time t) and the complementary function <j) c {t), defined as: 

Mt) = 1 - m (io) 

In the appendix we show that one can write the following equation for <p(t): 

m = /%c(o • p ( i - -r%T )dc (ii) 







MO 



m) M t) ■ P ( 1 - ] (12) 



that is: 



dt " w V 4>c{t), 

The failure time r is defined as the time at which all the fibers are broken. 
In the appendix we show that in the limit <f> « 1 and 2a/2 KT < (1 — f ) it is 
possible to write an approximated solution of eq.[l2| with a normal distribution 
of the thermal noise. This solution is: 



and 



^^(cy^, 



0(t)~-J,ln(l-i) (14) 

where 

c _ (1 - fo)fo 

KT ' 

Notice that, except for logarithmic corrections, the main dependence of r on 
fo and of KT in cq.H and in eq.O is the same. 
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3.2.2 The lifetime of the disordered fiber bundle 



In the disordered case the breakdown threshold f c is a random variable. In order 
to find an analytical solution of t as a function of /, one has to make several 
approximations, which are described in the appendix. When f c is normally 
distributed with mean < f c >= 1 and variance KTd = one finds that 

in the limit of </» << 1, 2V2 KT < (1 - f Q ) and 2^2 KT d < (1 - f Q ) a good 
approximation of the failure time is: 



with 



2 KT eff 



KT 

KT eff * — — (16) 



and 



2(l-/o) 



2 V2ttKT 



(17) 



One notices that for KTd ^ there is an effective temperature, which is an 
increasing function of the disorder variance KTd- This important observation 
will be numerically checked in sect. 4. 

Another interpretation of eq.|l5| and eq.|l^ is that the disorder changes the 
critical force from f c = 1 (in the ordered case) to a smaller value (1 — Y / 7rcr /2) in 
the disordered case. We consider that the interpretation in terms of KT e f f is the 
best one. Indeed one can show (see [pi}) that if the dynamics would be dominated 
by the fibers with thresholds in the tails of G G (/ C ) than KT eff = KT + KT d . 
The interpretation given in eq.^| stresses the difference between a tail dominated 
dynamics and a dynamics dominated by the fibers with f c > f m in- 

4 Numerical results 

In this section we check the analytical results obtained in sect. 3 on a numerical 
simulation on the fiber bundle. In the following we will assume without loss 
of generality f c = 1 and N = 1000. We have checked that the results do not 
change considerably up to N — 10 6 . For each values of the parameters KT, F, 
KTd we have repeated the experiments at least 10 times in order to estimate 
the scattering of the results in different realizations. 
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4.1 Failure time at constant force 



4.1.1 The homogeneous case 

In order to study the behaviour of failure time r, we started by imposing a force 
F(t) = F constant in time (creep test). If disorder noise and thermal noise 
are both zero (i.e. KT = KTd — 0), the system either breaks at time r = 
only if F > Nf c or it does not break at all otherwise. Conversely, if thermal 
noise is not zero (i.e. KT ^ 0), we observe (fig. || a) that failure time r has an 
exponential dependence on (1 — / ) 2 for any fixed value of KT, specifically: 

r ~ exp [a(l - f f] (18) 

where a is a fitting parameter, which is a function of KT. Furthermore at 
constant f a (see fig. || b), the failure time r depends on KT as follows: 

T ~ r ° eXp (^) (19) 

where A is a function of /. We notice that these results are similar to those 
predicted by the activation model discussed in section 3 and that the lifetime 
is well described by eq.||. More precisely one can check the predictions of cq.fL3| 
by plotting in fig.||c log(-^==) as a function of ■ ^ e see tnat a ^ data 

collapse on a single straight line of slope 1, verifying in this way the hypotheses 
made to obtain eq|l3[ In fig.|| we plot the time evolution of $ as a function of 
time. The continuous line is the numerical solution of cq.[l2| The circles corre- 
spond to the mean values obtained in several direct numerical simulations of the 
DFBM. In the inset we compare the result of the direct numerical simulation 
with the approximated solution (eq.jWJ) of eq.|l2|. 

Fig.|| and fig.|| show that the hypotheses made to obtain eq. H and eqjH are 
well verified by the numerical integration of the model. 



4.1.2 The heterogeneous case 

We are now interested in studying the dependence of the failure time r on 
disorder noise. To this aim, we have done simulations keeping KTd fixed at 
constant (non-zero) values. In fig.f| we plot the dependence of t / \Z2nKT, at 
fo = 0.54, as a function of 1/KT for different values of KTd- We notice that for 
the same values of KT, r decreases by several order of magnitude by increasing 
the disorder noise. Most importantly, we see that for any values of KTd the 
leading dependence of r on KT is the same of the ordered case (eq.|l9j) , but the 
value of A and r depend on KTd- Specifically, increasing KTd the value of A 
decreases whereas r Q increases. Using eq.|l^ and eq.[l6| one obtains that A in the 
heterogeneous case is: 

A = °- 5(1 -« 2 ( 1 -^)) 2 ™ 
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From the best fits of the numerical data of fig.^| one can measure the values of A 
as a function of / and KTd- The prediction of eq|2C| can be accurately checked 
by plotting 2 4/(1- f ) 2 as a function of ^2-rtKTd/ (1 — fo)- hi fig|5| we clearly 
see that the numerical values of 24/(1 — / ) 2 , obtained for different values of 
fo, collapse on the same curve, which agrees quite well with the theoretical 
prediction, continuous curve, obtained using eq|2C . 
The dependence of r on KTd is plotted in fig. 



for several values of / and 



KT. The continuous lines are the predictions of eq. 15 



We see that the main 



dependence of r on KTd is very well described by eq.|15|. The deviation at large 
KTd are due to the approximations made to get the analytical expression cq.[l5|. 
This equation and numerical results show that t q increases as a function of KTd. 
Therefore the observed reduction of the sample lifetime for increasing KTd is 
mainly due to a decrease of A. In looking at fig.f| and from eqjl^ we also notice 
that the derivative of t with respect to KT decreases as a function of KTd- 
Summarizing, the following conclusions can be extracted from the numerical 
and the analytical results: 

Bl The failure time r decreases as the disorder noise variance KTd increases, 
that is the more the medium is heterogeneous, the smaller the failure time 
is, see fig|| 

B2 As the disorder noise variance KTd increases, the derivative of r with 
respect to KT decreases, that is failure time r becomes less sensitive to 
the effective value of the thermal noise variance, as shown in fig.|]. 

Thus the role of the spatial disorder is that of multiplying KT by a factor, 
which is an increasing function of KTd- This is equivalent to say, as we already 
mentioned in section 3, that an approximation of KT e ff is given by eq.[l6|. 
Furthermore fig.|] and fig.|] show that the disorder reduces the dependence of r 
on the variation of the temperature T. 



4.2 Comparison with experimental results 

These results allow us to give a reasonable explanation of the experimental 
observations reported in the introduction. In fact, experiments show that the 
lifetime r of very heterogenous materials is independent of T fl6[] , while the 
lifetime t of quite homogenous materials as 2D-crystals depends on temperature 
rfll4|. This is consistent with the numerical results of our model. As shown in 
figJjT when KTd is small as in 2D-crystals the dependence of r on KT remains 
important, whereas this dependence becomes negligible when KTd is huge (like 
in chip-board panel wood and fiberglass). In any case disorder noise induces an 
effective temperature much larger than the thermodynamic one T. 



4.3 Failure time distribution 

We have studied the distribution of failure time r when the parameters KT and 
F are constants, in the case of a single fiber and in the case of a net of many 
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fibers. Clearly, in the case of a single fiber the probability of failure is constant 
in time, so that r follows a Poisson's law fig]?]. Conversely, in the case of a net 
of fibers one sees that r is normally distributed. This result is coherent with 
experimental results, which show that in highly disordered materials and large 
aspect ratio the sample lifetime distribution is Gaussian. 



5 Generalization to a time dependent imposed 
force 

In ref. ||l6f we proposed to generalize eqj^ such that it may be applied to any 
time dependent pressure; a similar procedure is proposed here. Suppose that 
is the damage density per unit time of the medium subjected to a force F . 
For an arbitrary load F(t) the total damage density at time t is: 



r* 1 

1(f) = / _ exp 

Jo T 



(l-F(t')/Nf 



eff 



dt' (21) 



thus the certitude of failure is attained when 1(f) = 1. Notice that this yields 
eqj| if F(t) = F. 

In order to verify whether this equation holds, we have done simulations impos- 
ing F(t) with different dependencies on time. 

As an example we consider a linear dependence on time of F: F — M t. 
Inserting this linear dependence in eq.^l] we can approximately compute the 

dependence of r on the slope M getting for M << 1 1 ,v 1 ■ 



t = 



1 - V^KT 5 



M [1 + 0.5 V2KT I 5 M 

where 

5 M -- 



(22) 




In fig.|| we plot r as a function of M at KT = 0.01. Circles are the computed 
values and the continuous line is the prediction of eq. |2^. The agreement 
between theoretical prediction and numerical simulations is quite good. 

Further we apply to the bundle a sinusoidal force as shown in fig.^. In this 
figure the dotted line is the external applied force and the dashed line repre- 
sents the numerical integration of the integral in eq.|2l]. The circles indicate the 
occurrence of breaking events. Finally the dotted line is the fraction of broken 
fiber <fi as measured by the direct numerical simulation of the fiber bundle. We 
clearly see that the lifetime predicted by eq.|l] perfectly agrees with that of the 
direct numerical simulation. This shows that eq.|2l] gives a good estimate of the 
breaking time for a periodic forcing. 
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6 Kaiser effect 



In 1950 Kaiser discovered that the acoustic emission of a stressed metal sample 
is zero if the applied stress is smaller than the previously applied maximum P8fl . 
This effect, usually called Kaiser effect, was also discovered in rock materialsp9|, 
|30f , but its existence was seriously questioned for the Westerley granite J3lj and 
for heterogeneous materials as wood and fiberglass Jll| . In our simulation we 
call event the simultaneous breaking of several fibers, and size of the event the 
number s(t) of fibers which break. The energy e associated to an event is the 
sum of the energies released by the fibers that crack, that is : 



s(t) ( F 



2Y \N-n{t) 



(23) 



The potential energy of a fiber that breaks is proportional to its accumulated 
elastic potential energy. If we assume that the acoustic energy is proportional to 
this potential energy then e may be compared to the acoustic energy measured 
in the experiments [fl0[ |ll|| . The cumulative energy E a (t) is the sum of e from 
to t. 

In order to check the behaviour of our model with respect to the Kaiser effect 
we have imposed cyclic forces F = Fq[1 — cos(^ L i)] to the net (the continuous 
line in fig. |l0|a). If KT = 0, i.e. in a classic DFBM, the Kaiser effect is strictly 
valid: no events appear after the first period. Conversely, if KT > 0, events 
(spots on fig. [Io|si) are produced after the first. In fig. [To|b), the number of 
events Ni (continuous bars) and the energy Ei (dashed bars) released during 
the i-th cycles are plotted, normalized to the values of the first cycle, as a 
function of the cycle number i. The released energy E2 and the number of 
events N2 occurred during the second cycle are not negligible as they represent 
about the 50% of E\ and N\. From the second cycle (i > 2), Ni and Ei grow 
with i. During the fourth cycle the energy released £4 becomes bigger than 
the energy idreleased during the first cycle. The number of events 2Vj grows 
slower than the released energy Ei so that the average energy Ei/N released 
by an avalanche during the i-th cycle grows as a function of i. The behaviour 
of Ei /N is represented by circles on fig. |ic|b) . The rates of the 2Vj and Ei after 
the first cycle depend on KT and KTd- The statistic of events is very similar 
to experimental observations in heterogeneous materials JHJ. In conclusion, we 
have shown that the Kaiser effect is valid only if the thermal noise is negligible. 
Conversely, if a thermal noise is present the Kaiser effect is not observed. 



7 Conclusions and discussion 

To simulate the failure of heterogeneous materials, we used a DFBM with ther- 
mal noise. We observe that the breakdown of the system occurs through the 
coalescence of many microcracks which are driven by the thermal fluctuations 
and which appear at different times. Thus the lifetime r of the DFBM with 
thermal noise follows a law, which is similar to an activation processes. The 
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main goal of the simulation was to check the role of the disorder on the lifetime 
of the system. We have shown both numerically and analytically that the most 
important effect of the disorder is that of producing an effective temperature 
of the system, which is larger than that of the thermal noise. This effective 
temperature (eq.|l^) is equal to the thermal noise temperature multiplied by a 
factor, which is an increasing function of the disorder variance. This important 
result is independent on the specific distribution of thresholds. In this sense 
we can say that the disorder "amplifies" the thermal noise temperature. It is 
important to stress that that thermal noise and disorder play different roles: 
indeed it is the thermal noise, which drives the cracking process, while the dis- 
order just enhances its effect. Two important hypotheses have been done to 
get this result. The first is that the threshold distribution evolves as a function 
of the fraction of broken bonds. The second is that the tails of the distribu- 
tion do not play an important role in the fracture process. To estimate r, one 
has to consider only the part of the distribution containing about 90% of the 
unbroken bonds. The effective temperature of eqjlf] is quite different from the 
result recently obtained by Roux^I]. He finds that the effective temperature 
is the sum of KT and K Td- Indeed our result is not directly comparable with 
that of Roux because he estimates the time needed to cut the first fibers. This 
time is negligible with respect to the total lifetime of the bundle. Furthermore 
the number of bonds contained in the tails does not affect the resistance of the 
system (for example in the case of the Gaussian only 2% of the bonds have a 
threshold in the tails). 

Our results not only demonstrate that the disorder enhances the role of 
thermal noise but eqjll] and fig.^ clearly show that when the disorder noise 
variance KTd increases, the absolute value of the derivative d d £ T decreases, that 
is the failure time r becomes less sensitive to variation of the effective value of 
thermal noise (not zero as well) . The conclusion, that in an activation processes 
disorder noise induces an effective temperature larger than the thermodynamic 
one, has been reached for other disordered systems such as foams This 
kind of disorder induced enhancement of the thermodynamic temperature is an 
interesting result because it could be be a quite general property of disordered 
systems, where a thermal activated processes with long range correlation is 
present. 

It is interesting to compare the results of our model with the recent experi- 
mental observations on the lifetime of samples submitted to a constant applied 
stress. It has been shown that in wood fiberglass Jl6| , gels jl5| and 2D- 
crystals the dependence of the lifetime on the applied stress is well predicted 
by the Pomeau model on the delayed fracture (see eq.|l|). However we already 
mentioned in the introduction that in all of the experiments the estimated tem- 
perature was higher than the thermodynamic one and in ref. we argued that 
the heterogeneity of the material may enhance the role of temperature. The 
theoretical and numerical results of this paper give new insight to this argu- 
ment. The fiber bundle and the Pomeau model have a different dependence on 
the applied load (compare eq.|| and eqjl]). This difference comes from the fact 
that the energy barrier is different in a real crack and in the bundle. Thus we 
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neglect this difference and we focus on the role played by the disorder. The 
Pomeau model is based on the idea that failure of the sample is due to the ther- 
mal nucleation of one defect (microcrack) initially present inside the material 
[ p"7| . In contrast, this model and experiments show that the failure is due to the 
coalescence of many microfractures appeared at different times. We think that 
a good theoretical model to describe the dynamic of the DFBM and various 
real materials submitted to an external force could be made transforming in a 
statistical model the one of Pomeau. In that model the nucleation of a micro- 
crack is the elementary process, i.e. the sample break when a certain number 
of microcraks have nucleated. In this case, the parameters Y, T and T e // of 
eq]l], become average parameters. The strong time-dependent fluctuations of 
the internal forces induced by the heterogeneity (defects, microcraks ...), can 
be considered as a sort of noise so that, in our statistical model, T e ff must 
depend on the thermal temperature but also on the disorder in the medium. 
As we have already mentioned our numerical and analytical results on the fiber 
bundle also show that the dependence on the lifetime on temperature decreases 
when the disorder increases. These results could explain why in very disordered 
materials, such as wood and fiberglass, the lifetime r seems to be independent 
of temperature while r strongly depends on temperature in weak heterogeneous 
materials such as 2D-crystals. 

The comparison of the lifetime of the bundle with those of real materials 
merits a special comment. Indeed in the bundle we consider as unit of the time 
a complete update of all the fibers. One can assume that the characteristic 
time in a real material is the time needed to redistribute the stress after that a 
microfracture is occurred. This time, which is given in a first approximation by 
the size of the sample divided by the sound speed, is of the order of 10~ 5 s 
This means that the largest observed experimental time correspond to roughly 
10 10 time steps of the bundle. This very large r are reached for very small KT 
and KTd where our approximation becomes more precise. Thus we may assume 
that the comparison of this theoretical results with experiment could be rather 
realistic. 

In the last part of the paper we have introduced the concept of density of 
damage per unit time finding an equation (eq.^l]), which predicts the lifetime r 
of the DFBM for any time dependent imposed force. These numerical results 
agrees with experiments on wood and fiberglass Finally, we have shown 
that the validity of Kaiser effect is linked to the value of the thermal noise KT 
and of the disorder noise KTd- 
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APPENDIX 
Analytical results 



A The homogeneous case 

Consider the homogeneous case where the threshold force f c is the same for all 
of the N fibers composing the bundle. Let 



P( V ) 



f 00 f 00 1 / x 2 \ 



be the probability that the thermal noise has a fluctuation larger than 77. If n 
is the number of broken fibers at time t then the force present on each fiber is 
/ = /oj^j' Thus the probability that a single fiber breaks in a time step is: 

Pl = P(f c -f) = P (f e (A.2) 
The expected number of fibers that break in a time step is: 

R N = (N-n)-p 1 = (N-n)-p(f c - 1 ^Sj ( A .3) 

Notice that if this expression is smaller than one (i.e. i?jv < 1) one may think 
of (t) = as the expected time between two consecutive microcracks. 

These expressions allow us to get several analytical results. Indeed it can be 
easily proved that as N — > 00 the ratio becomes independent of N: 



i=l 

= s£<*-»^(/c-^)= (A.5) 



1=1 

t 



We can introduce the function <p(t) (fraction of broken fibers at time t) and the 
complementary function 4> c (t), defined as: 

Ht) = lim | (A.7) 



s (t) = 1 - m (a.8) 
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Turning the sum in eq. A. 6 into an integral one can write the following equation 
for (f>(t): 



(1 - 0(C)) ■ P [fa - 



fo 



i - HO 



that is: 



dcj>{t) 
dt 



(i - m) p[fc- 



fo 



i - 4>(t) 



(A.9) 



(A.IO) 



The failure time r is defined as the time at which all of the fibers are broken. 
It can be expressed in terms of <p(t) by the condition: 



1 



(A.ll) 



In the homogeneo us cas e we assume f c = 1 for all of the fibers. Unfortunately 
no exact solution of ( A.IO ) is known, so one has to make some approximations 
to get an analytical result. The crucial point is that in the direct numerical 
simulation of DFBM one observes that <j> <C 1 for most of the time: this allows us 
' and to approximate P(rf) using an asymptotic development: 



to write 



1- 



P(V) = o 1 - erf ( 



'\/2fcT' 



Thus one can write: 
d4 

~dt ~ 



kT 1 



*(i-A) 

kT \-<t> 
2^1-fo-fo 

being f Q (j) <C 1, one may write: 
t~kT 1 



exp 



■ exp 



IkT 1 



1-0 

2kT 



(1 - fo ~ lo<i 

IkT 



d0 
~dl 



2tt 1 - f 
= B exp(C0) 

being B and C defined as: 



exp 



(w r 

2kT 



exp 



2(1 - fo)U 
2kT 



B 



C = 



kT 1 
V 2^ 1 - fo 

(1 - fo)fo 

kT 



exp 



(l~/o) 2 

2kT 



(A.12) 
(A.13) 

(A.14) 
(A.15) 
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At this point one turns to the integral equation: 

J exp (-C<f>)d<j) = B J dt 
The last equation can be integrated: 



cxp(-C^) R , , 1 

— c — = - Bt+ c 



(A.16) 



where the initial condition is (f>(0) = 0. 



We verified numerically that this equation describes the system behaviour 
for small <f>. Moreover one finds that as <f> becomes big enough, say <f>break — 0.1, 
the system undergoes a cascade of mi crocracks until the final breakdown, which 
occurs after a little time. Eq.( A. 1(f ) seems to represent closely the system's 
behaviour for almost all the time. 



Thus one can write the following estimation for the failure time t: 



C B 



■ exp 



(Wo): 

2kT 



(A.17) 



This expression of r agrees with both numerical simulations and theoretical re- 
sults of H. 



Using ( |A.17|) into ( |A.16| ) one finds 

T - t 



exp (— C <j)) 



(A.18) 



Finally one finds: 



In 



r - t 



-C 4> 



(A.19) 



We verified (see fig.§) that these equations agree very well with numerical data, 
actually even for values of 4> comparable to 1. From eq.A.18 one obtains that 



t ~ 0.95r at <j> = 3/C. Thus the system undergoes a fast cascade of microcracks 
until the final break up for 4> > tftbreak = 3/C: 



3 KT 



^break 



3 (1 - fo) 



2irKT\ 



(A.20) 



fo(l-fo) 2 /„ (lnr-0.5 ln-^ 
This equation verifies that (j> << 1 for any value of f a if KT is small enough 
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B The disordered case 



In the disordered case the breakdown threshold f c of fibers is a random variable. 
Let G(f c ) be the distribution of thresholds, that is G(f c )df c is the fraction 
of fiber with threshold bigger than rj. The crucial point is that the distribution 
of f c evolves with time (i.e. (f>) as fibers break, in such a way that: 



G{f c A)df c = i-4> 



The expected fraction of fibers which break at a given time t is: 



(B.l) 



Rd 



P[f, 



fo 

1-c 



G(f c ,4>) df c 



We use the definition of P(r]) (see eq.A.l) to write: 



R<t, = / / G T (x) G(£,<j>) dt]dx 

poo pOO 

/ GrUc-OG^^d^dU 

fo 



(B.2) 



(B.3) 



where f' a - 

The function if)(f c , = Gr{f c — f ) G(£, <j>) is integrated in the plane (/ c , ^) 
within the region Re which is defined by the intervals £ < f c < oo and < £ < 



co. In cq.(B.3) one may invert the order of the integration by considering that 
the inte rvals defining region Re arc: < £ < f c and < f c < oo. Therefore 
cq.(B.3) becomes: 



Rd> — 



Gt( fr - 



fo 



dfc 



fc 



Thus in the disordered case one may turn eq.( A. 10 ) into 
Gt (f, 



dt 



(B.4) 



(B.5) 



with the condition imposed by eq.B.l. Notice that if the a delta distribution 



with normalization B.l is chosen for G then one gets from cq. B.5 the equation 
for the homogeneous case eq.(A.lO). Equation eq. B.5 is far too complex to 
be integrated analytically. However one can greatly simplify this framework 
by assuming that the dynamics is controlled by the weakest fibers which are 
the first to break. This hypothesis can been checked in a direct numerical 
simulation using for Go(f c ) a Gaussian and an uniform distribution centered 
at f c = 1. The result of the numerical simulation are shown in fig.|ll| for the 
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Gaussian distribution. Wc clearly see that are the weakest fibers which break 
first and that the distribution evolves just on the left side. 

Starting from this numerical evidence of the behaviour of G(f c ) one can find 



an approximated analytical solution of eq.B.5 



B.l Case of the normal distribution for thresholds 

Let the initial distribution Go(/ c ) of thresholds be a normal distribution of 
variance kTd — o%/2 and mean value (/ c ) = 1. In fig{n] we see that ( for 
2V2KT < (1 - f ) ) when increases the distribution is deformed on the left 
side, whereas the right side remains unchanged. On a first approximation one 
can take into account this behaviour by splitting the distribution at time t: 



G(fc) 



■ exp 



(fc 



for f c > 1 



(B.6) 



and 



G(f c ) = exp ( - (/c 2 ly I for f c < 1 (B.7) 



J TV (To 

where Od depends on (f>. The normalization condition |B.1| yields ad — a a (1 



2(f>) , which roughly agrees with the numerical data. Thus the integral in eq. B.4 
becomes 



1 

^ jo 
1 

7t 



exp 



(fc l-ch) 2 



1 - erf( 



1 - fc x \ &d df, 



<*d 



2 a a 



exp 



(fc 1 °j > ) 



' °d . t Jc - 1 A df c 
— +crf( ) — 

O n On AO 



(B.C 



In order to find an analytical expression for eq. B.8 one has to consider that in 



the case of a Gaussian distribution 90% of the unbroken fibers have a threshold 
fc > fmin = (1 — a d/2). One can make the hypothesis that only these 
fibers, with a large threshold, play an important role in supporting the applied 
load. The mathematical counter part of this hypothesis is that only the part of 
the distribution with f c > f m in h as to be taken into account. Thus in the first 
integral of eq.B.8 the lower integration limit becomes f m in = (1 — \[k ^d/2) = 
[1 — ^/tt <t (1/2 — <p)] and we can use erf(x) ~ 2x/y/ir . After applying these 
approximations one gets: 



R * = ^ L„ exp [ — — ) l^w^ ) — + 

1 r ( (fc-^? \ df c 

-= exp (1-0) — B.9 
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Taking into account that, within the integration limits of the first integral, 



(fc — fmin) is almost constant with respect to the variation of exp I — 



one can take the mean value of f c — /„ 
approximated result: 



Therefore one obtains the following 



1 



erf 



So 



erf 



(B.IO) 



In this equation one can use the asymptotic approximation of the error 
function erf(x) ~ 1 — ex^p(—x' 2 )/(x^/^^). Thus one gets the equation describing 
the time evolution of d>. 



dt 



B n (l-$exp (A n (f>) 
where the coefficients A n and B n arc the following: 

(1_ v^£o -/ )(/ -V5F<t 



A, 



kT 



(B.ll) 



(B.12) 



B " = kT 



exp 



2^ (l + exp(- ^<!- /} ) 



2(1- 



fo) 



fo 



2kT 



(B.13) 



Following the same procedure used in the homogeneous case, one finds that 
the failure time is: 



2 v / 2^fcT 



(fo - a )(l + exp(- ^f^Q )) 
"(l-/o) 2 



exp 



2 feT 



1 - 



/7T CT n 



2 (1 ~ /o) 



(B.14) 



Notice that for ad — > one recover eq. A. 17 for the life time in the homoge- 
neous case. In section (4) we have seen that the r predicted by eq.( B.14 ) agrees 
very well with the results of the numerical simulations. 



B.2 The uniform distribution 



The same calculation can be done taking for G(f c ) in eq. B.l a uniform distri- 
bution of mean value < f c >— 1 and width D. We find in this case that: 



21 



VMT 2 

(f - 0.8 *£>)(! + exp(-^^Sj) X 



exp 



(Wo) 2 ^_ 0.8 D \ 2 ' 



2 kT V 2(1- /„) 



(B.15) 



where the factor 0.8D comes from the fact that in the case of the uniform 
distribution 90% of the the unbroken bonds have f c > / m i n (l — 0.8 D/2). 

C Discussion 

One observes that in both cases of uniform and normal distributions, the pres- 
ence of a disorder on thresholds have the same effect over the general behaviour 
of the system. In particular one sees that disorder enhances the effect of ther- 



mal noise, in that it decreases in a multiplicative way the exponents in B.14 and 



B.15. We stress that thermal noise and disorder play different roles: indeed it 



is the thermal noise which drives the cracking process, while the disorder just 



enhances its effect. The effective temperature of eq.16 is quite different from the 
result recently obtained by Roux. He finds that the effective temperature is the 
sum of KT and KTd- Indeed our result is not directly comparable with that of 
Roux because he estimates the time needed to cut the first fibers. To compute 
this time Roux considers that the threshold distribution is almost quenched. 
This is certainly correct at the very beginning of the dynamics but it is not 
at longer times as we have shown in previous sections. In any case the time 
needed to cut the first fibers is negligible with respect to the total life time of 
the bundle if (1 — f Q ) > \J2 KTd- Furthermore the number of bonds which have 
the threshold in the tails of the distribution does not affect the resistance of 
the system, for example in the case of the Gaussian only 2% of the bonds have 
the threshold in the tail and in the case of a uniform distribution there are no 
tails. 
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b) 




Figure 1: a) Modified DFBM: N fiber in parallel, with the edges fixed on a 
rigid support are subjected to an externally imposed force F which is distributed 
democratically on the net, that is all fibers (not broken) are affected in the same 
way. Each fiber is also subjected to a random (zero mean, normally distributed), 
addictive force A/j(i), where is intended that A/i(i) is a realization of a white, 
time independent stochastic process. Name n(t) the number of broken fibers 
at time i; we derive the following expression of local force fi for the i-th fiber: 
fi = N F n ^ + A/i. b) The equivalent of the fiber bundle model is the fuse 
network. One can think to our model as a fuse network where the Nyquist noise 
of several resistances (current generators 5Ii) are the noise generators for each 
bond 
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(a) (b) 




5 10 15 



(1-f//(2 KT) 

Figure 2: Homogeneous net (KTd = 0) in a creep test, (a) Failure time r 
as a function of (1 — / D ) 2 for several values of the thermal noise variance KT 
((x),KT = 0.0045, {o),KT = 0.006, (v),ifT = 0.01). The continuous lines 
are the fits with eq.[l8| (b) r as a function of 1/ KT for several values of imposed 
force f-N. (□)/« = 0.45; (o) f a = 0.54; (v) f a = 0.7. The continuous lines are 
the fits with eq.||. (c) r' = r / (2tt KT)- 1 / 2 as function of (1 - f ) 2 /(2KT). 
The same symbols of (b) are used. The continuous line is obtained from eq.[13|. 



24 



10 1 10 2 10 3 1 4 

Time 

Figure 3: Time evolution of cj> at / = 0.6, KT = 0.008 and KT d = 0. The 
continuous line corresponds to a solution obtained by the numerical solution of 
cqjl^. The dashed line is instead the approximated solution (eq|l4|) of eq.p^|. 
The symbols (o) correspond to the results of the direct numerical simulation of 
the DFBM. These data are plotted in the inset as a function of 1 — t/r in a 
semilog scale. 




40 60 80 100 120 140 160 180 200 



1/KT 

Figure 4: Heterogeneous bundle (KTd ^ 0). Dependence of the lifetime r 
on 1/KT at f — 0.45. The different symbols correspond to different values of 
KT d : (o) KT d = 0, (<) KT d = 0.02, (D)KT d = 0.04 Notice that both r and 
J^ T decrease as KTd increases, that is the more media become heterogeneous, 
the smaller is r and the dependence on KT of r. 
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Figure 5: Dependence of 2.4/(1 — f ) 2 (A is defined in equation |2^) as a 
function of y/2ir KT d /(l - f a ) for different values of f a : (□) f a = 0.45; (o) f = 
0.54; (*) f Q — 0.7. The continuous line is the theoretical prediction obtained 
from eq.|2(]. The symbols (v) correspond to a numerical simulation done with 
uniform distribution at f Q — 0.54. In this case the x axis is 0.8D/(1 — f )( see 
appendix B.2). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 6: The lifetime r is plotted as a function of the disorder noise vari- 
ance KTd for several values of KT and of f a and different threshold distri- 
butions. The continuo us lin es represent predictions of eq.[l5] for the Gaus- 
sian distribution and eq. B.15| for the uniform distribution. Discrete points are 
the results of numerical measures. For the Gaussian distribution: (o) KT = 
0.01, f a = 0.54; (□) KT = 0.01, f a = 0.45; (o) KT = 0.007, f Q = 0.54; (< 
) KT = 0.01, f a = 0.6; (*) KT = 0.0073, f Q = 0.45. For the uniform distri- 
bution (+) KT = 0.01, f = 0.54; (y) KT = 0.01, f Q = 0.45. 
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Figure 7: Lifetime distribution for a single fiber (a) and (b) for N 
So = 0.54 and KT = 0.021 in both cases 



1000. 




Figure 8: Failure time r as a function of the slope M in the case of a linear (in 
time) imposed force F — M t. The theoretical prediction of eqp2](continuous 
line) confirms the results obtained in numerical simulations. In these simula- 
tions: N = 1000, KTd = 0, (o) KT = 0.005 and (□) KT = 0.02 
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Figure 9: A sinusoidal force (dotted line) is applied to the bundle. The dashed 
line represents the numerical integration of the integral in eq.^l]. The circles 
indicate the occurrence of breaking events. Finally the continuous line is the 
fraction of broken fiber <p as measured by the direct numerical simulation of 
the fiber bundle. The time predicted by eq.^l] perfectly agrees with the that 
of the direct numerical simulation. In this simulation N = 1000, KTd = 0, 
KT = 0.003 
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Figure 10: (a) Imposed force F = F a [l — cos(^-t)} as a function of time in a 
single realization; discrete points represent the occurrence of breaking events, 
(b) Number of events N{ (continuous bars) and the energy E 1 , (dashed bars) 
released during the i-th cycles, as functions of i. Both these quantities are nor- 
malized to the values of the first cycle. Discrete points represent the normalized 
mean energy of an event occurred during the i-th cycle. We observe that 

Kaiser effect is not valid, that is there are a large number of events after the 
first cycle. The parameters of this simulation are N — 1000 and KT = 9 ■ 10~ 3 
KT d = 1.5 • 10~ 2 . 
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Figure 11: Heterogeneous net. Evolution as a function of <f> of the threshold 
distribution at / = 0.6,/TT = 0.006 and KTd = 0.005. (a) Evolution of G(/ c ) 
at (ft = (o),(ft = 0.05 (o),(f> = 0.1 (+) and (ft = 0.15 (□). The fraction 
4> = 0.15 corresponds to t ~ 0.9r. (b) Evolution of the integral distribution 
^(/ < /c) = Jo" G(x)dx at the same (ft used in (a). 
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Answer to referee A: 

We thank the referee for his useful comments. The manuscript has been 
changed in order to keep them into account. Specifically the answers to the 
questions of the referee are the followings: 

0) We agree the abstract contained too strong statements. It has been 
changed. 

1) All these suggestions have been taking into account 

2) We thank the referee for having indicated these references which were 
unknown to us. References to the work of Coleman and Phoenix have been 
added. 

3) We agree the statement has been changed. 

4) We agree. The validity limits of our approximation are now clearly 
indicated in the text. 

5) We added the symbol for / = 0.7 in the figure. 

6) These changes have been done. 

7) We agree. The text has been corrected. 

8) The suggestion has been taken into account. 

9) We agree is the fraction of fiber. It has been changed. 

10) Now we explain how to go from eq.B2 to eq.B3 (now eq.B4). 

11) We agree with the referee that our assumptions were badly explained 
in the previous manuscript. This part has been completely rewritten. We 
agree that our approximations do not work for large temperature and this is 
now clearly indicated. However we disagree with the referee on two points. 
The fact that the threshold distribution is not quenched and that the fibers 
with threshold larger than one are not broken are not assumptions but they 
come from a numerical evidence. We use this evidence in order to find an 
analytical formula which works for the cases we have analyzed. Furthermore 
the 90% has nothing to do with temperature because it is reasonable to think 
that the most important fibers are those with the largest thresholds. Indeed 
this approximations works quite well also for the uniform distribution which 
has no tails. From this point of view there is no reason to believe that our 
approximation fails for T going to zero. However we completely agree with 
the referee on another point. The comparison with the results of Roux's 
paper was incorrect and imprecise. The discussions in appendix C and in the 
conclusions of the paper have been changed. 
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Answer to referee B: 

We thank the referee for the careful reading of the manuscript. His sug- 
gestions have been taken into account and a figure with cycling loading has 
been added. The paper has been carefully read. 

The comments of the referee allows us to make our point of view very 
clear. 

The recent experimental results were the motivation of our work. We 
wanted to study the influence of disorder in the crack nucleation processes. 
The fiber bundle is a certainly a toy model which could be useful to under- 
stand what is going on. To make this point clear we changed the abstract 
and part of the introduction. 

However we strongly disagree with the referee when he says that every- 
thing is clear and the only unexplained experiments are ours. We kindly 
explain the reasons: a) The temperature is very high in all the experiments. 

In our experiment it is 3000-fT In the Pauchard experiment (ref.14) the tem- 
perature is larger than 1000 K. In the Bonn experiment (ref.15) it is larger 
than 10 W K. The referee can easily check our claims taking the data of ref.14 
(PRL of Pauchard et al.) and ref. 15 (Science, Bonn et al.) In order to make 
this point clear we wrote these temperatures in the introduction. 

b) The value of the dimensionality d, measured in Bonn's paper, is ob- 
tained from fig. 4 of the paper on Science (ref.15). In this figure the vertical 
axis is the log of the log of the lifetime. However by making the second log 
the authors of ref.15 forgot that there is an additive constant (the log(r )). 
If the calculation is done taking into account the log(r ) (this can be approx- 
imately obtained using the values given in the text of ref.15) one gets (if a 
power law still exists) values of d which are very different from those claimed 
in ref.15. Thus there is something wrong in that paper. 

c) The Pauchard data are much better fitted by the Mogi's law than the 
Pomeau's law. A simple chi-square test shows this. Further I want to recall 
to the referee that in that paper only the best three points have been chosen. 
"Best" means here the points which give the best fit on the Pomeau's law 
and not the points with the smallest statistical errors. This kind of data 
analysis can be found in a paper published on PRL. 

d) Poissonian statistics are observed only by Bonn. In the other experi- 
ment the life time distributions are Gaussians. 

We did not add points b) and c) in the introduction of our paper because 
it will generate an useless and nasty polemics with Bonn and Pauchard. We 
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would like to avoid this polemic. However if the referee thinks that it is 
important to put the points b) and c) in the introduction of our paper, we 
are ready to do that. We think that this is not really a key point and it will 
make our introduction too long. We want just to kindly drive the attention 
of the referee on the fact that in spite of his nasty comments our experiment 
is the best verification of the Pomeau's law, as far as the force dependence 
is concerned. At this point the referee has two alternatives. Either he think 
that Pomeau's law is useless or that it is an extremely good idea but one 
has to improve the understanding of its applicability. This is exactly what 
we think and what we have done in our paper. We show that disorder can 
enhance the role of thermal noise at least in a toy model. Of course the 
activation energy of the fiber bundle is not the same of a continuous medium 
and we discuss this point in section 3.1. However the thermal activation of 
the fiber bundle seems to be already so reach which certainly merits to be 
studied. 
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